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abstract. In 8 , Rugang Ye proved the existence of a family of constant mean 
curvature hypersurfaces in an m + 1-dimensional Riemannian manifold (M m+1 ,g), 
which concentrate at a point po (which is required to be a nondegenerate critical 
point of the scalar curvature), moreover he proved that this family constitute a 
foliation of a neighborhood of pq. In this paper we extend this result to the other 
curvatures (the r-th mean curvature for 1 < r < m). 
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1 Introduction 

Let S be an oriented embedded (or possibly immersed) hypersurface in a Rieman- 
nian manifold (M n+1 , g). The shape operator As is the symmetric endomorphism 
of the tangent bundle of S associated with the second fundamental form of S 1 , 65, 

by 

b s (X,Y)=g s (A s X,Y), VI,FeT5; here g s = g\ TS . 

The eigenvalues Ki of the shape operator As are the principal curvatures of the 
hypersurface S. The fc-curvature of S is define to be the k-th symmetric function 
of the principal curvatures of S, i.e. 



H k (S) := ^2 ^ 



■ <ik 



Hence, when k = 1, Hi is equal to n times the mean curvature of S. When 
M n+1 = W l+1 is the Euclidean space, H2 is equal to n ("~ 1 ) times the scalar cur- 
vature of S and H n is equal to the Gauss-Kronecker curvature of S. In this paper 
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we are interested in the existence of hypersurfaces in M n+1 whose fc-curvature is 
constant. Hypersurfaces with constant mean curvature, constant scalar curvature 
or constant Gauss-Kronecker curvature in Euclidean space or space forms constitute 
an important class of submanifolds. In Riemannian manifolds very few examples of 
constant /c-curvature hypersurfaces are known, except when k = 1. 

R. Ye H3, has proved the existence of a local foliation by constant mean 
curvature hypersurfaces which concentrate at a point (which is required to be a 
nondegenerate critical point of the scalar curvature function). We extend the result 
and methods of [5| to handle the case k = 2, . . . , n. No extra curvature hypotheses 
are required. In particular, we prove the existence of foliations of a neighborhood 
of any nondegenerate critical point of the scalar curvature of (M n+1 ,g) by constant 
Gauss-Kronecker or constant scalar curvature hypersurfaces. As in the idea is 
to perturb S p (p), a geodesic sphere with small radius p > centered at a point p. 
A simple computation will show that S p (p) is close to being a constant k-curvature 
hypersurface as p tends to and in fact 

a k (S p (p)) = C*p- k + 0(p- k + 2 ), 

In this paper, we show that it is possible to perturb S p (p) for every small radius, 
to a constant /c-curvature hypersurface equal to C k p~ k for any 1 < k < n — 1, 
provided p is close to a nondegenerate critical point of the scalar curvature of M. 
The analysis here is inspired from the one performed in (Hj. In fact, independently of 
the value of k, the linearized /c-curvature operator about the unit Euclidean sphere 
is always a multiple of Agn + n, the linearized mean curvature operator about the 
unit Euclidean sphere. This implies that, as in to perform the perturbation 
of a small geodesic sphere, one has to overcome the problem of the existence of 
(n + l)-dimensional kernel of Ag™ + n, kernel which is related to the invariance of 
/c-curvature with respect to the action of isometries (in the case of the unit sphere, 
this kernel is only generated by translations). This is where, as in [5] we use the 
fact that we are close to a nondegenerate critical point of the scalar curvature of 
the ambient manifold. 

We notice that the analysis performed in |S] is specific to treat the case of mean 
curvature, namely k = 1 and, unfortunately, can't be used to treat the general case 
k = 2, . . . , n. The main technical result of this paper is a precise expansion of geo- 
metric operators (first and second fundamental forms) for perturbed geodesic sphere 
(see Proposition 12. II Proposition 13. II and Proposition 13. 2[) . We believe that these 
expansions are of independent interest and can be used in many other construction 
0]. Our main result is : 

Theorem 1.1. Suppose that pg is a nondegenerate critical point of the scalar cur- 
vature 1Z of M . Then there exists po > 0, such that for all p e (0, po), the geodesic 
sphere S p (po) may be perturbed to a constant k-curvature hypersurface S p with 
Hk = C k p~ k . Moreover these k-curvature hypersurfaces constitute a local folia- 
tion of a neighborhood of pa . 
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The existence of the hypersurfaces is not so difficult and can be obtained rather 
easily. The fact that they constitute a local foliation requires more work. The leaves 
S p are small perturbation of geodesic spheres in the sense that S p is a normal graph 
over Sp(po) for some function w p which is bounded by a constant times p 2 . 

The hypersurface S p is a small perturbation of S p (po) in the sense that it is 
the normal graph of some function (with L°° norm bounded by a constant times 
p 2 ) over a geodesic sphere obtained centered at a point at distance bounded by a 
constant times p 2 of pq. 

Existence of families of constant mean curvature hypersurfaces concentrating 
along positive dimensional limit sets is obtained by R. Mazzeo and F. Pacard in j^] 
and then in collaboration with the author in [2] in a more general setting. 

The paper is organized in the following way: In section 2 we expand the coeffi- 
cients of the metric in normal geodesic coordinates. Section 3 will be devoted to the 
expansion of the first fondamental form, second fondamental form and the Shape 
operator of the perturbed geodesic spheres. Using these, we derive in section 4, the 
expansion of the fc-curvature of the perturbed spheres. Section 5 is devoted to the 
proof of the main result of this paper, theorem ll.il 

2 Expansion of the metric in geodesic normal co- 
ordinates 

In this Section we introduce geodesic normal coordinates in a neighborhood of a 
point p E M . We choose an orthonormal basis E{, i = 1, . . . , n + 1, of T P M . 

Consider, in a neighborhood of p in M, normal geodesic coordinates 
F(x) := ex.pp 4 (xi EJ, x := (xi, . . .,x n+ i), 

where exp M is the exponential map on M and summation over repeated indices is 
understood. This yields the coordinate vector fields Xi := F*(d Xi ). As usual, the 
Fermi coordinates above are defined so that the metric coefficients 

9ij = 9(Xi,Xj) 

equal Sij at p. We now compute higher terms in the Taylor expansions of the 
functions <7y. The metric coefficients at q := F(x) are given in terms of geometric 
data at p := F(0) and |x| := (x\ + . . . + ^ + i) 1/2 - 

Notation The symbol C(|x| r ) indicates an analytic function such that it and its 
partial derivatives of any order, with respect to the vector fields x- 7 Xi , are bounded 
by a constant times \x\ r in some fixed neighborhood of 0. 

We now give the well known expansion for the metric in normal coordinates 1 , 
|7|, but we will briefly recall the proof in the Appendix for completeness. 
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Proposition 2.1. At the point q = F(x), the following expansions hold 

9ij = Sij + ^g(R(E k ,E i )E e ,E j )x k x e (2.1) 

+ \ g(V Ek R(Et, Ei)E m ,Ej)x k x e x m + 0{\x\ 4 ) 

where all curvature terms are evaluated at p. Convention ever repeated indices is 
understood. 

3 Geometry of spheres 

In this Section, we derive expansions as p tends to for the metric, second funda- 
mental form and mean curvature of the sphere S p (p) and their perturbations. 

Fix p > 0. We use a local parametrization z — > Q(z) of S n C T p M. Now define 
the map 

G(z) :=F(p(l-w(z))Q(z)), 

and denote its image by S p (p, w), so in particular S p (p, 0) = S p (p). Because of the 
definition of these hypersurfaces using the exponential map, various vector fields we 
shall use may be regarded either as fields along S p (p, w) or as vectors of T p M. To 
help allay this confusion, we write 

Q — ^Ej Q t := d z iO j Ej. 

These are all vectors in the tangent space T p M. On the other hand, the vectors 

T — QiXj Ti:=d z iQ j Xj 
lie in the tangent space T q M , where q — F{z). For brevity, we also write 

Wj := d z jw, Wij := d z i d z jw. 

In terms of all this notation, the tangent space to S p (w) at any point is spanned 
by the vectors 

Z j = G*(d z j)=p((l-w)r j -w j r), j = l,...,n. (3.2) 

Notation for error terms 

The formulas for the various geometric quantities of S p (p,w) are potentially very 
complicated, and so it is important to condense notation as much as possible. Fortu- 
nately, we do not need to know the full structure of all of these quantities. Because 
it is so fundamental, we have isolated the notational conventions we shall use in 
this separate subsection. 

Any expression of the form U (w) denotes a linear combination of the functions 
w together with its derivatives with respect to the vector fields 0j up to order j. 



4 



The coefficients are assumed to be smooth functions on S n which are bounded by 
a constant independent of p £ (0, 1) and p £ M, in C°° topology. 

Similarly, any expression of the form Q J (w) denotes a nonlinear operator in the 
functions w together with its derivatives with respect to the vector fields Oj up 
to order j. Again, the coefficients of the Taylor expansion of the corresponding 
differential operator are smooth functions on S n which are bounded by a constant 
independent of p £ (0, 1) and p £ M in the C°° topology. In addition Qi vanishes 
quadratically at w = 0. 

Finally, any term of the form LP x Q k will denote any finite sum of the product 
of a linear operators L J with nonlinear operators Q k . 

We also agree that any term denoted 0(p d ) is a smooth function on S n which 
is bounded by a constant (independent of p) times p d in the C°° topology. 

The first fundamental form The next step is the computation of the 
coefficients of the first fundamental form of S p (p,w). We set q := G(z) and p := 
G(0). We obtain directly from QTfy that 

giX^Xj) = 5 lJ + \g{R{Q,E l )Q.E j )p 2 {l-w) 2 

+ \g{V & B{Q 1 E i )Q,E j )^0.-wf (3.3) 

+ O(p 4 ) + p i L (w)+p 4 Q°{w). 

where all the curvature terms are evaluated at p. Observe that we have 

0(T,T) = 1 ffpr.T^sO 

Using these expansions it is easy to obtain the expansion of the first fundamental 
form of S p (p, w). 

Proposition 3.1. We have 

p- 2 (1 - w)- 2 g(Zi, Z : ) = 3(6,, 9,) + | g(R(e, Q t ) 6, 8,) p 2 (1 - w) 2 

+ i 3 (v i?(e, eo e, e,-) P 3 (1 - w) 3 + (1 - w)- 2 w t Wj + o{ P A ) 

+ p i L°(w)+p i Q°(w). 

(3-4) 

where all curvature terms are evaluated at p. 

The normal vector field Our next task is to understand the dependence 
on w of the unit normal N to S p (w). Define the vector field 

N := - y + A> Zj, 

and choose the coefficients so that TV is orthogonal to all of the Zj. This leads 
to a linear system for A 1 . 

^ A J g(Zj, Zi) = -puii 

3 
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Observe that 

g(N, N) = l + pJ2 A i w j 

3 

The unit normal vector field N about S p (p, w) is defined to be 

N 

N := (3.5) 

g(N,N)V* 

The second fundamental form We now compute the second funda- 
mental form. To simplify the computations below, we henceforth assume that, at 
the point e(z) 6 S n , 

g(Qi,Qj) = 6ij and Ve.Oj = 0, i,j = l,...,n (3.6) 

(where V is the connection on TS 1 " -1 ). 

Proposition 3.2. The following expansions hold 

-g(V Zt N,Z 3 ) = p(l-w)S ij +pw ij + y(R(e,e i )Q,Q j )p 3 (l-w) 3 

+ ^.g(v e i?(e,e l )e,e,)p 4 (i-^ 4 

- | (g(R(Vw, Qi) 9, 6,) + g(R(Q, 9,) Vw, 6,)) p 3 
+ 0( P 5 ) + p A L\w) + pQHw) + P L 2 {w) kQ\w) 

(3.7) 

where as usual, all curvature terms are computed at the point p. 

Proof : We will first obtain the expansion of g(\7 Z .N, Z 3 ). To this aim, we 
compute 

= g(W Zi T, Zj) g{V Zi (A k Z k ), Z 3 ) 

= t^5(Vz, ((1 - w)T), Z 3 ) + T ^ Wi g(T, Z 3 ) <?(Vz 4 (A" Z k ), Z 3 ) 
= T^5(Vz, ((1 - w) T), Z 3 ) -j^w t Wj 9(V Zi (A k Z k ), Z 3 ) 

Now, recall that 

Y / A k g{Z k ,Z j ) = -pw 3 
k 

Hence 

Y,9(VzM k Z k ),Z 3 ) = -pw t3 -Y J A k g(Z k ,V Zt Z j ) 

k k 

Using the fact that 

2g(Z k ,'Vz i Zj) = Z i g{Z k ,Z J ) + Z 3 g(Z k , Zi) - Z k g(Zi,Z 3 ) 
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we conclude that 

9(Vz< (A k Z k ),Z 3 ) = -p Wij ~\Y, Ak (Zi9(Zk, Zj) + Z j9 (Z k , Zi) 



2 

k k 

- Z k g{Z l ,Z :j ) 

To analyze the term V^((l — w) T), let us revert for the moment and regard w 
as functions of the coordinates z and also consider p as a variable instead of just a 
parameter. Thus we consider 

F(p, z) = F(p(l-w(z))e(z)). 

The coordinate vector fields Zj are still equal to F*(d Zj ), but now we also have 
Zq := (1 — w) T = F*(d p ), which is the identity we wish to use below. Now, we 
write 

g(V Zi ((l - w) T), + g(V Z] ((1 - w) T), Z«) 
= g(V Zi Z ,Z j )+g(V Zj Z ,Z i ) 
= Z g(Zi,Zj) 

Collecting the above we have obtained to formula 

-ffCVzjiV, Z,-) = ji^) Z a 9( z h Zj) - j^pWi Wj + pwij 

+ I E fc ^ 9(Zk, Zj) + Zj g(Z k , Zi) - Z k g(Z h Zj)) 
We will now expand the first and last term in this expression. 

If the coordinates y are chosen so that g(Oi,Qj) — 5™ at the point where we 
will compute the shape form, we have, using the result of Proposition l3.il 

^^ZogiZ^Zj) = p(l-w)5 lJ + lg(R(e,O t )e,e 3 )p 3 (l~w) 3 

+ £ <?(v e i?(e, e,) e, e 3 ) P i (1 - w f + T ^ pWl Wj 

+ O(p 5 )+p 5 L (w)+p 5 Q°(w). 

Using the same Proposition together with the fact that the coordinates y are chosen 
so that VQ.Qj = at the point where we will compute the shape form, we also 
have 

Zi g(Z k , Zj) + Zj g(Z k , Zi) - Z k g(Zi, Zj) = 
f (g(R{Q k ,Q i )e j ,e)+g{R(Q k ,e i )Q j ,e)) p 4 

+0(p 5 ) + p 2 L\w) + p 2 Q 1 {w) + p 2 L 2 {w) k L\w) 

If the coordinates y are chosen so that g(Qi,@j) = 8ij at the point where we will 
compute the shape form, we have the expansion 

A " = v +pL 1 (w) + P Q 1 (w) 

p(l — wy 
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collecting the above estimates, we conclude that 

-gptztNiZj) = P (l-w) Sa +p Wij + | g(R(9, 8,) 6, 8,) p 3 (1 - w) 3 

+ ^g(v e R(9,9 i )9,9 j ) P i (i-w) 4 

- I {g{R{9 k , 8,) 8, 8,) + g{R(9, 8 4 ) 9 k , 8,)) p 3 w k 
+ 0( P 5 ) + p A L x (w) + pQ\w) + p L 2 (w) x Q 1 (w) 
It remains to observe that 

g(N,N)- 1 / 2 = l + Q 1 (w) 

This finishes the proof of the estimate. □ 

The shape operator of perturbed surfaces Collecting the esti- 
mates of the last subsection we obtain the expansion of the shape operator of the 
hypersurface S p (p,w). In the coordinate system defined in the previous sections, 
we get 

Proposition 3.3. Under the previous hypothesis, the shape operator of the hyper- 
surface S p (p,w) is given by 

pA l0 {w) = (1 + w)Sij +Wij + %g(R(Q, 9i)9,9j)p 2 

- | [g{R{9, 8,) 8, 8,) w + g(R(0, 0<), 8, 9 k ) w kj 

+ (g(R(9 k , fy) 8, 8,) + g(R(9, 8 4 ) 9 k , 8,)) w k ] p 2 

+ \ 5 (V e i?(8, 80 8, 8^) p 3 + G(p 4 ) + p 3 L 2 (w) 

+ Q 1 (w) + L 2 (w) k L°(w) + L 2 (w) k Q l {w). 

where all curvature terms are computed at the point p. 

4 The fc-curvature of the perturbed sphere 

Given any symmetric matrix A, and any k = 0, . . . , n, we define 

°k{A) := K---K- 

ix<...<i k 

where Ai, . . . , A„ are the eigenvalues of A. The k-th Newton transform of A is 
defined by 

T k (A) := a k (A) I - a k . x (A) A + -- + (-l) k A k . 

with T n (A) = 0. Now suppose that A — A(t) depends smoothly on a parameter t, 
it is proved in [a] that 

j t a k {A)=Tt(r k _ 1 {A)j t A} (4.8) 



8 



From this computation, it follows at once that, given any n x n symmetric matrix 
H, 

a k (I + H) = C k n + G k -_\ Tr(H) + 0(\H\ 2 ) 

Using this together with the previous expansion of the shape operator, it is not hard 
to check that the fc-curvature of the hypersurface S p (p, to) can be expanded as 

p k H k (S p (p,w)) 

= Ct + C k n Z\ [(A S » +n)w-\ Ric(9, 6) p 2 - \ V e Ric(9, 9) p 3 
+ | (Ric(9, 9) + 2 Ric(V-, 9) - g(R(@, V-)9, V-)) w p 2 

+ 0(p 4 ) + p 3 L 2 (w) + Q x {w) + L 2 (w) k L°(w) + L 2 (w) K Q^w)] 
where as usual, all curvature terms are computed at p. Here we have defined 

Ric(V-,9) := Ric(e l ,9)e l 

and 

g(R(0, V-), 9, V-) := g(R(Q, e 4 ), 9, e,) e 4 e 3 

if ei, . . . , e„ is an orthonormal frame field of TqS n satisfying V ei ej = at the point 
q G S n where these expressions are computed. It will be convenient to set 

C := 1 (Ric(9, 9) + 2 Ric(V, 9) - g(R(Q, V), 9, V)) 

Now observe that a similar expansion is valid in Euclidean space and in this case 
the expansion of p~ k H k {p,w) does not depend on p (nor on p). This means that 
the nonlinear operator 

Q 2 := Q 1 + L 2 k L° + L 2 x Q 1 

can be decomposed into its value in Euclidean space and a similar operator all 
of whose coefficients are bounded by p. This fact can also be recovered by going 
through all the above expansions. Therefore, we can write 

Q 2 = Ql + P Ql 

where Q 2 is the corresponding nonlinear operator when the metric is Euclidean and 
hence it does not depend on p; while p Q 2 denotes the discrepancy induced by the 
curvature of the metric g on M. Both Q 2 and Q 2 satisfy the usual properties. 

5 Existence of foliations by constant /c-curvature 
hypersurfaces 

Assume that we are given p S M, a nondegenerate critical point of the scalar 
curvature 1Z on M. We would like to find a small function w G C 2 ' a (S n ) and a 
point p close to po such that 

H k (S p (p,w))=C k p- k 
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In view of the previous expansion, this amount to solve the nonlinear equation 

(Agn + n) w = I Ric(e, 6) p 2 + ~ V e Ric(6, 6) p 3 - 0{p 4 ) 

6 4 (5.9) 

-p 2 £w-p 3 L 2 (w)-Q 2 (w) 

We denote by II and II the L 2 -orthogonal projections of L 2 (S n ) onto Ker(Ag>i + 
n) and Ker(Agn + n) , respectively. Recall that the kernel of Agn + n is spanned 
by (fi, for i = 1, . . . , n + 1, the restriction to the unit sphere of Xj, the coordinates 
functions in R" +1 . 

First fixed point argument From now on, we assume that the function w G 
C 2 ' a (S n ) is L 2 -orthogonal to Ker(Agn + n) and we project the equation (|5.9(l over 
Ker(Ag™ + n) ± . We obtain 

(A g n + n) w = n x [| Ric(e, e) p 2 + i v Ric(e, e) P 3 - o(p 4 ) 

-p 2 £w-p 3 L 2 (w)-Q 2 H] 
We define wq & Ker(A5>. + n) 1 - to be the unique solution of 

(A S n +n) Wo = iRic(e,9) (5.10) 

since 11^ (Ric(9, 9)) = Ric(9, 9). Similarly, we define wi £ Ker(A s « + n) x to be 
the unique solution of 

{A S n + n) Wl = i n x [V e Ric(9, 9)] 

It is easy to rephrase the solvability of the nonlinear equation l|5.9|l as a fixed 
point problem since the operator Ag™ + n is invertible from the space of C 2 ' a (S n ) 
functions which are L 2 -orthogonal to Ker(Ag« + n) into the space of C°' a (S n ) 
functions which are L 2 -orthogonal to Ker(Agii + n). We write w := p 2 wq + p 3 w\ + 
p v, so that it remains to solve an equation which can be written for short as 

(A S n + n) v = -0(1) - p- 2 Cw - p- 1 L 2 (w) - p- 4 Q 2 (w) 

Applying a standard fixed point theorem for contraction mappings, it is easy to 
check that there exists a constant n > 0, which is independent of the choice of 
the point p £ M , such that there exists a unique fixed point in ball of radius n in 
C 2,a (S n ), provided p is chosen small enough, say p G (0,po)- We denote by v p this 
solution and define 

w p :— p 2 Wo + p 3 w x + p 4 v p . 
It is easy to check that, reducing the value of po if this is necessary, 

\\wp~Wp'\\ C 2. a{Sn) < cp 2 dist(p,p'), (5.11) 
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for some constant c which does not depend on p £ (0, po) nor on p or p' . In addition, 
the mapping 

(p,p) G (0,po) x M — > ttfp G C 2 < Q (S") 

is smooth and 

||-DpWj,|| C 2, Q(s „ ) +p||5pWp|| C 2. Q (5>.) < cp 2 

for some constant c which does not depend on p G (0, po) n or on p. 

Second fixed point argument It now remains to project the equation (|5.9|) 
where w has been replaced by w p , over Ker(As» + n). To this aim, we recall the 
nice and key observation from [S]- 

The problem is to compute the L 2 -projection of the term <?(Ve-R(6, Oi)0, Qj) 
over the kernel of the operator A$n + n. This amounts to compute, for any m — 
1, . . . , n + 1, the quantity 

Xj Xk Xl x m 



Now to evaluate this quantity, simply use the fact that the integral vanishes unless 
all indices are all equal or constitute two pairs of equal indices. Using this, together 
with the symmetries of the curvature tensor which imply that R(E, E) = 0, we 
obtain 

B m = g{VE m R{Ei,E m )Ei,E m ) (J Sn xf-3 J s „xlx 2 2 ) 

+ g{V Em R(E i ,E j )E i + 2V Ej R{E i ,E m )E i ,E j ) f s „ x\ x\ 
Now, use second Bianchi identity 

g{W Em R{E l ,E J )E l ,E J ) = 2g{V Ej R(E i ,E m )E i ,E j ) 
together with the fact that 

/ 4_o / 2 2 _ ^ f 2 

/ x \ ~ 6 I x l x 2 ~ 7 T~^T / x l 

is- Js*> (n + 3) J Sn 

To conclude that 

n ( 5 (v e i?(e, Qi) e, e,)) = — l —g{vn, Xi e { ) 

n + 3 

where 1Z denotes the scalar curvature function, computed at p. 

Therefore, the projection of the equation l|5.9|l over Ker (Ag™ + n) yields 

g(Vn,x i E i )=V p 

where we have defined 

V p := 4 (n + 3) II [p~ 3 0(p 4 ) + p' 1 C w p + L 2 {w p ) + p~ 3 Q 2 (w p )] 



11 



Now, using the fact that po is a nondegenerate critical point of the scalar curvature, 
we conclude easily (applying for example a topological degree argument) that there 
exists p close to pq satisfying (|5.11|) provided p is close enough to 0. This gives 
the existence of constant fc-curvature leaves for all p small enough, unfortunately it 
turns out that the point p is at most at distance a constant times p from po and 
this is not enough to show that the constant fc-curvature leaves form a foliation of 
a neighborhood of pq. 

To improve this estimate, many observations are due. First, observe that we 
can decompose 0(p 4 ) into the sum of two functions, one of which is homogeneous 
of degree 4 (in the coordinate functions Xi) and the other one which is bounded 
by a constant times p 5 . The L 2 -projection of the homogeneous function of degree 
4 is equal to since this homogeneous function is invariant under the change of 
coordinates 6 into — 0. Hence we conclude that 

|n(eV))|<cp 5 

Similarly, observe that Wq and hence C Wq are invariant under the change into 
— and hence the L 2 projection of Cwo over Ker(As»i +n) again identically equal 
to 0. Therefore, we conclude that 

\U(Cw p )\ < cp 3 

Finally, we use the observation at the end of §4. Since the nonlinear operator Q 2 
preserves functions which are invariant under the action of —J, we conclude that 
U(Ql(p 2 wq)) = and hence 

|n(Q 2 K))| <c P 5 

These precise estimates imply that, 

\V p \<cp 2 

for some constant which does not depend on p nor on p. With slightly more work, 
we get using similar arguments that 

|II (y p - V p >) I < cp 2 dist(p,p') (5.12) 

Now, for all p small enough, we can find a solution of (|5.9|l using a fixed point 
argument for contraction mapping, in the geodesic ball of radius 2 p 2 centered at 
any nondegenerate critical point of 1Z. Moreover, the solution p p depends smoothly 
on p and 

\d P Pp\ <cp 

This later fact, together with (|5.12(l shows that the solutions constitute a local 
foliation. This completes the proof of the main result. □ 

Having derived such precise estimates, we can compute the expansion of the n- 
dimensional volume of the leaves of the foliation as well as the (n + l)-dimensional 
volume enclosed by each leaf. 
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Proposition 5.1. For all p small enough the following expansions hold for the 
n- dimensional volume of S p 

Vol n (S p ) = p n Voln(S n ) (1 - i^y^p 2 + o( P 4 ) 

and the (n+ 1) -dimensional volume of the set B p enclosed by S p and containing the 
point po 

1 / 77 4- 2 

Vol n+1 (B p ) = — — p n+1 Vol n (S n ) 1 - —L—TZp 2 + 0(p 4 ) 

77+1 \ 277(77 + 6) 

where the scalar curvature is computed at Pq, a nondegenerate critical point oflZ. 
Proof : Integrating (|5.1U|I over S n we find 



J n w ° = \ ^„ Ric ( e ' 9 ) 



Now, plugging the expansion of w p into the expression of the first fundamental form 
given in Proposition 13. II we find the expansion of h the induced metric on S p 

p- 2 h l3 = (1 - 2p 2 w - 2p 3 Wl )S l3 + 1 g(R(Q, Q z ) 9, 8,) p 2 

(5.13) 

+ ±g(V B R(e,Q i )e,Q J )p 3 + 0(p 4 ) 

This implies that 

P~ n V\h\ = l-np 2 w -np 3 w 1 -^ Ric(9, 6) p 2 - V e Ric(9, 9) p 3 + 0(p 4 ) 

The first estimate follows from integrating this expansion using the fact that the 
integral of wi and the integral VeRic (9, 9) over S n vanish together with the fact 
that 



/ Ric(9, 9) = — !— Vol n (5 n ) n 

J S" 77+1 



Next, we consider polar geodesic normal coordinates (r, 9) centered at p p . In 
these coordinates the metric g expanded as 

r~ 2 gij = Sij + J g(R(Q, 9,) 9, 9,) r 2 + 1 5 (V e i?(9, 9,) 9, 9,) r 3 + 0(r 4 ). (5.14) 

then, the volume form can be expanded as 

r- n V\g\ = 1 - I Ric(9, 9) r 2 - i V Ric(9, 9) r 3 + 0(r 4 ) 
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Integration over the set r < p (1 — w p ) give 

voi„ +1 (£g 



r <P (1— Wp) 



r n I 1 - ~ Ric(6, 9) r 2 - 4: V Ric(6, 6) r 3 + 0(r 4 ' 
6 12 



' -P n+1 I (l-(n + l)p 2 w ) + O(p n+5 ) 



n+1 Js 

lJ_ p n + S 

6 n + 3 J s 



(1 - (n + 3)p 2 w )Ric(e,6) 



1 r 1 n™+ 3 , 

-p' l+1 Vol n (5") - p n+3 / w - - ^— / Ric(6, 0) + 0( P n+5 ) 



n+1 J S n 6 n + 3 j g 

_L^»+iVol n (S") - ( 1. + 1 — j— )p»+ 3 / Ric(0, 9) + 0(p™+ 5 ) 
n + 1 3n b n + 3 Jon 



1 ■ VoIn(S") f 1 - + 0(P 4 ) 



n + 1 r v ' V 2n(n + 3) 

This gives the second estimate. This proves the desired result. □ 

6 Appendix : proof of Proposition 12.11 

The aim of this Section is to prove Proposition 12. II Observe first that the curve 
s — > expp f (si?) is a geodesic. Therefore, if X is the unit tangent vector to the curve 
we have V xX = 0. Hence we also have (V x) m X — for all m > 1. In particular, 
we have, at p, 

(v E y n E = o 

for all E G T p M and for all m > 1. 

Observe that X a are coordinate vector fields hence 

Vx a Xb = V ' x b X a 

Taking E = E a + e E b and looking for the coefficient of e in V E E — 0, we get 

V Ea E b = 

Looking at the coefficient of s in V 2 E E — 0, we get 

2V 2 Ea E b + V Eb V Ea E a = (6.15) 
Finally, looking at the coefficient of e in V^-E = 0, we get 

2 V 3 E b + (V Bi V £a + V Ea V Eb )V Ea E a = (6.16) 
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Recall that, by definition 

V x V Y := R(X, Y) +V y V x + V: x ,y; 
Hence, if X and Y are coordinate vector fields we simply have 

V x V Y X := R(X, Y)X + Vy V X X (6.17) 

We also have 

VyV^VyX := Vyi?(x,y)x + i?(VyX,y)x + i?(x, v Y y)x 

(6.18) 

+ R(X,Y)V Y X + V 2 Y V X X + VyV [x , y] X 
Now use (|6.15|l and i|6.17|) to obtain 

3V 2 Ea E b = R(E a ,E b )E a , (6.19) 
Similarly, use l|6.16|l and l|6.18|l to obtain 

2V Ea E b + R(E b ,E a ) V Ea E a + 2V E y Eb V Ea E a =0 (6.20) 
Since V Ea E b = 0, we get 

2 V| a £ h + 2V Ba V Ei V Ea E a = 
Using this, we conclude that 

2V 3 Ea E b = -2V Ea V Eb V Ea E a = -2V Ea {R{E b ,E a )E a + V Ea V Ea E b ) 

= -2V Ea {R{E b ,E a )E a )-2V Ea E b 

Hence 

2 V| a S fc = -V Ea R{E bl E a ) E a (6.21) 

Now, we have 

X c g ab = g{V x a X a ,X b ) + g{X a , V Xc X b ), 
and we get X c g ab \ p = 0. This yields the first order Taylor expansion 

g ab = <U + 0(M 2 ), 

To compute the second order terms, it suffices to compute X 2 g ab at p and polarize. 
We compute 

X'lgab = g(V Xc X a ,X b )+g(Xa,V Xc X b ) + 2g(V Xc X a ,V Xc X b ) 
Using (|6.2U|) we get 

XZg ab \ p = ^g(R(E c ,E a )E c ,E b ). 
15 



The formula for the second order Taylor coefficient for g ab now follows at once. 
Similarly, we compute 

X?g ab \ p = g(V Xc X a ,X b ) + 3g(W Xc X a ,W Xc X b ) 
+ 3g(\7 Xc X a , V 2 x X b ) + g{X a , V x X b ) 

and using 1)6. 21[) this gives 

X*g ab \ p = g(\7 Ec R(E a ,E c )E b ,E c ). 
the formula for the second order Taylor expansion for g ab holds at once. □ 
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